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1 Introduction 



Mukhanov and Chibisov ^Lj were the first to suggest that quantum fluctua- 
tions during inflation produced the tiny inhomogeneities needed to form the 
various cosmic structures we observe currently - as the result of gravitational 
collapse over the course of more than 10 billion years. Early work on the sub- 
ject was also done by Hawking |2., by Guth and Pi [S], and by Starobinskii 
jl]. The formalism has since been described at length in a number of review 
articles OIHIIZI. It has received much attention recently owing to the unprece- 
dented precision with which the imprint of these fluctuations on the cosmic 
microwave background radiation has been imaged by the WMAP satellite 

Much of the fascinating structure revealed by these measurements derives 
from processes which occurred long after the end of inflation, and are not the 
subject of this paper. Instead, we re-compute the primordial fluctuation spec- 
trum which is the starting point for the analysis of subsequent processes. The 
justification is that we now have at our disposal the exact scalar and graviton 
mode functions upon which the calculation is based fOlE]- There has never 
been any doubt regarding the spacetime dependence of the mode functions 
during the epoch of matter domination in which the cosmic microwave back- 
ground radiation anisotropies accumulate. What was previously unavailable 
is an exact expression for the normalization factor which the mode functions 
build up during inflation. 

Previous computations have been based on approximation schemes that 
were developed over the course of two decades. A key step in this effort 
was the introduction, by Stewart and Lyth, of the slow-roll Bessel function 
approximation jT2]. However, Wang, Mukhanov and Steinhardt |TSj demon- 
strated that carrying this approximation to higher orders does not generally 
improve accuracy, while Martin and Schwarz ^3] showed that the technique's 
accuracy is not sufficient for comparison with precision experiments such as 
WMAP and PLANCK. Recent improvements [W\ [T7| ITH] have overcome 
these obstacles, at least for slow-roll inflation so the additional precision 
available from our exact solutions is probably not necessary for comparison 
with foreseeable data. But it is nice to have, and it is simple enough to con- 
struct exotic models in which the slow-roll paradigm breaks down completely. 
We shall study one in an appendix. 

To fix notation, note that cosmologically relevant spacetimes are charac- 
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terized by scale factor a: 



ds^ = — dt"^ + a^{t)dx-dx . (1) 

Although not directly an observable, the ratio of its current value Oq to its 
value at past time t is the cosmological redshift experienced by light emitted 
at that time and received now: 

m . - 1 . (2) 

Its logarithmic derivative defines the Hubble parameter H which measures 
the rate at which distant matter is receding due to the expansion of the 
universe: 

Its second time derivative enters into the deceleration parameter q: 

^^^^ - d^{t) ~ H\t) • 

The weak energy condition implies that q{t) > —1; inflation is characterized 
by q{t) < 0. 

Quantum fluctuations are not especially big during inflation, but they 
are enormously larger than afterwards. Therefore, we can analyze the pro- 
cess using linearized quantum field theory. Furthermore, the high degree of 
homogeneity and isotropy of the inflationary geometry implies both that a 
fluctuation can be characterized by its constant, co-moving wave vector k : 

S = ^ . (6) 

and that each fluctuation evolves independently. The physical wave length 
Aph of a fluctuation grows as the universe expands: 

Aph(^) = a{t) X , (6) 

where A is the co-moving wave length. During inflation the Hubble radius 
r-H- 

rnit) ^ H-\t) , (7) 
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Figure 1: The first and second horizon crossings of the physical wavelength Aph- The 

Hubble radius th is constant during inflation (blue), behaves like th ^ 
during radiation (red)^ and like th ~ a 2 during matter domination (green). The present 
is at oq. The graph is not properly scaled. 

is approximately constant, whereas it grows more rapidly than the scale 
factor after the end of inflation (see Figure 1). This variation of rn gives rise 
to the two horizon crossings which characterize the fluctuations of interest 
to us. They happen when: 

Horizon Crossing =^ \b{t) = fuit) . (8) 

First horizon crossing occurs during inflation. Before this time the linearized 
fields oscillate with falling amplitude; afterwards they are approximately con- 
stant. Second horizon crossing occurs long after the end of inflation, indeed 
after the emission of the cosmic microwave radiation. Before second horizon 
crossing the fields are approximately constant whereas they oscillate with 
falling amplitude afterwards. 

We can be more precise by defining the dimensionless variable x which 
represents the physical wave number in Hubble units: 



and in terms of which horizon crossing means: 



Horizon Crossing 



x{t, k) 



1 



(10) 



When the deceleration parameter q is constant, which is a good approxima- 
tion during the dominant phases in the history of the universe, the following 
relation is valid: 



If we take the initial time instant ti to signify the onset of inflation, it be- 
comes apparent that during inflation x decreases with time: q is negative in 
this era while z is an ever decreasing function of time. Thus, the modes of 
interest start with x larger than 1 and achieve condition (jlUp - first horizon 
crossing - as the universe still inflates. After first horizon crossing, the vari- 
able X for these modes further decreases and becomes very much less than 
1. However, post-inflationary evolution is characterized by positive q and, 
hence, x increases during this era. The modes of physical interest are such 
that condition ()10|1 is satisfied again - second horizon crossing - before the 
present. 

Moreover, we note that significant fluctuations do not occur for all fields, 
only for those which are both much lighter than the inflationary Hubble 
parameter and also not conformally invariant. These two requirements mean 
we need consider only gravitons and light, minimally coupled scalars. 

It is unnecessary to discuss invariant characterizations of cosmological 
perturbations. The fully general and invariant formula of Sachs and Wolfe 
- which is reviewed in Section 2 - allows us to solve for the perturbations 
with any convenient choice of gauge and field variables. Although we shall 
not work beyond linearized order, it is worth noting that the result of Sachs 
and Wolfe can be extended to any desired order in the weak field expansion. 
The method is applied for the generic system of a graviton with a massless, 
minimally coupled scalar in Section 3. The scalar and tensor power spectra 
are derived in Sections 4 and 5 respectively. In both cases improved esti- 
mates are obtained. Our conclusions comprise Section 6. The basics of the 
evolution dependent improvement factors have been summarized in the Ap- 
pendix. Another Appendix describes a model in which the slow roll paradigm 
completely breaks down but our methods can still be employed. 




(11) 
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2 The Sachs and Wolfe Effect 



The gravitational field equations are: ^ 

G^y = — -g^ivR — SttG T^i, , (12) 

in which G is the Newton constant. A spatially homogeneous and isotropic 
universe can be conveniently represented by the stress tensor T^^ of a perfect 
fluid with energy density p, pressure p and 4- velocity u^: 

Tf^u = {p + p)u^Ui, + pg^u , (13) 

where u'^ obeys: 

u'^u'^g^, = -1 . (14) 

To account for the observed structures and obtain a more realistic cosmo- 
logical description, deviations from homogeneity and isotropy are essential. 
Without any reference to their origin, it is simple to incorporate such depar- 
tures as linear perturbations on the dynamical variables of the system: 

gnv{r],x) = gt.ui'n) + 5g^^{ri,x) , (15) 
p{rj,x) = p{r}) + 5p{rj,x) , (16) 
p{r],x) = p{r]) + 5p{r],x) . (17) 

The unperturbed metric field g^^, belongs to the Robertson- Walker class of 
spacetimes and is, therefore, conformally flat and characterized by scale fac- 
tor a: 

g,jiv = a^{v)Vnu ■ (18) 

The unperturbed p and p correspond to the average energy density and pres- 
sure of the physical system respectively. The arbitrariness in the choice of 
coordinates is resolved by employing a frame that moves with the fluid: 

goo(v,x) = -1 , (19) 
G',irj,x) = ^ u^irj,x) = a-\rj) 5'', . (20) 

It is important to note the relation between co-moving and conformal time 
intervals: 

dt = a{rj) dri . (21) 
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Figure 2: The light emission {rjE , x) and reception (77^ , 0) events, and the hghthke 
geodesic x'^- 



Sachs and Wolfe have computed "511] the redshift accumulated by a light 
ray as it travels in the presence of (fT3j) from its emission to its reception (see 
Figure 2). The result is quite general as the only relevant ingredient is the 
metric field perturbation: 

Sg^,u = a^iv) h^y{r],x) . (22) 

If the light signal is observed from direction e, the wavelength shift z{e) is 
given by: 

' + ^ wka; • '''' 

where E and R stand for the emission and reception events respectively, and 
where k'^ is the 4-momentum of the light ray: 

k^ir) = rir) . (24) 



^By Rfj,i, and R we denote the Ricci tensor and Ricci scalar constructed from the 
spacelike metric tensor 17^^. Furthermore, an overdot indicates differentiation with respect 
to co-moving time t while an overprime denotes differentiation with respect to conformal 
time 77. 
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The lightlike geodesic x'^ satisfies: 

r{r) + n^[x(r)] r{r) X^{r) = , (25) 
9aMr)] r{r) x^{t) = , (26) 

and ()25II26|) can be integrated to give the following result for z{e) to first 
order in the perturbation h^y 



1 + z(e) 



X 



(27) 



1 - 



Vr~Ve 



da 



^■n _ ^ riB,~o, ere ) 



Suppose that thermal radiation of average temperature Te was emitted 
from a spacelike surface at the time tie of the coordinate system. Then, at 
the reception event (?7r,0): 



TR{e) 



TEie) 



(28) 



so that the first order temperature fiuctuation observed from direction e is 



ATe 
Te 

dt 



+ 



(29) 



tE 



dt' 



i(t') 



and has been expressed in terms of the co-moving time t. Since the accu- 
mulated wavelength shift and the temperature fiuctuation are observables, 
expressions fl27p and ()29|) are manifestly gauge invariant. 

The Cosmic Microwave Background Radiation (CMBR) consists of pho- 
tons emitted during the period of decoupling. ^ What is measured is the 
product of temperature fiuctuations simultaneously observed from two dif- 
ferent directions (see Figure 3). Thus, the connection between the measured 
quantity and the quantum mechanical origin of these fiuctuations comes from 
the study of: 



for the appropriate vacuum state \ fl). 



162 J 



n 



(30) 



^In the history of the universe, the period of decoupHng is centered around z ~ 1089 
with a width Az ^ 195. 
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observer 



• • • - z 

1200 1000 

Figure 3: A typical lightlikc geodesic on its way to the observer. The graph is not 
properly scaled. 

3 The Perturbations of the Gravit on- Scalar 
System 

The Lagrangian describing the system of the graviton and a minimally cou- 
pled scalar is: 

^ = T^R^a - Id.cpd^^g'^^V^ - V{v)V^ . (31) 

iOTrCr z 

Its dynamical variables are the metric field g^^i: and the scalar ip. Both are 
expressed background plus a quantum field: 

9t,u{v,^) = a^{v){v,iu + ' (32) 

= a^{r]) (^r]f,^ + n^^^{ri,x)^ , (33) 
ip{ri,x) = ifioiv) + 0(^,^) , (34) 

where k"^ = IGttG is the loop counting parameter of quantum gravity. 
The background Einstein equations are: 

3H' = 8nG{^^l + V{vo)) , (35) 
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{-l + 2q)H' = 8nG{^ipl-VM) , (36) 

where an overdot represents differentiation with respect to the co-moving 
time t. ^ Ahhough it is traditional to regard the potential as known and 
then infer the scale factor, with our method it is more convenient to regard 
a{t) - and hence H{t) = a/a and q{t) = —ad /a? - as a known function from 
which the background scalar and the potential can be expressed as follows: 



^« = -\/rT^ ^ . (37) 

V(»>„) = (2 - q(t)) ^ . (38) 
We parameterize the third derivative of a{t) using the variable: 



Hence, the derivative of the potential is: 



V\^,) = ( 2 - q[t) + r(t) ) + q{t) ^ . (40) 

Higher derivatives of the potential can obviously be obtained by taking higher 
time derivatives of the scale factor, for example: 

^"(^o) = ( 3g(t) r{t) -r\t)-^^ H\t) . (41) 



A convenient diagonalization of the linearized system is given in [211 122] 
and is summarized in pjj. By employing a generalized de Bonder gauge 
condition: 



F^, = a 



2 

all linearized fields can be expressed in terms of: 



, (42) 



V^5^(r/,x) = 5: {e,,(^,s)^^(r^,fc) e^^^-a(fc,.) + (c.c.)} , (43) 

^oo(r/,x) = j-^, [Uc{r^,k)e'^'-'Yik) + (c.c.) } , (44) 



■^The relation between co- moving and conformal time derivatives is: ^ = ^; 



as follows: 



ipoii'n,x) = , 



(f){r],x) 



1 1 d 



fYTm Hit) ait) dt 
The mode functions Ua,c are of the form: 



— (a(t) tpooiv,^) 



(45) 
(46) 

(47) 



ait) 



QAiv^k) , 



UAiv.k) = 



Uciv,k) = -^l + q{t)Hit)-Qciv,k) 



(48) 
(49) 



where Qa,c obey [TT] : 



Q 



A,C 



+ 



Oa = a 



1C 



Qa,c = 

1 



a vT+g 



(50) 

(51) 



The graviton and scalar creation and annihilation operators are canonically 
normalized: 



aik, s) , a\k', s') 
Yik) , Y\k') 



= i2n)' 5\k - k') 5ss' , (52) 
= (27r)3 6^ik - k') . (53) 

The graviton polarization tensor is purely spatial, transverse and traceless: 

eo^ik,s) = kieijik,s) = euik.s) = . (54) 

Moreover, summing products of two polarization tensors gives: 

1 



J2^ijik,s) e*^^ik,s) 



Iljj — 5ij ki kj , 



(55) 
(56) 



^ While the general solutions to (|50I51|) are known 11 , it is only in a particular limit 
that we shall need them for the purposes of this paper. 
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where Hij is the transverse projector. 

In order to apply the basic formula (|29|) for the first order temperature 
fluctuation, we must transform the hnearized fields ()44II47|) from obeying the 
gauge condition (jl^ to satisfying the co-moving gauge conditions P9II20|) . 
This is achieved by effecting the field-dependent coordinate transformation: 



which imposes the co-moving gauge conditions: 



eo{r],x) 
ei{r],x) 



dr]' a{r]') ^ooiv'^x) , 



dr]' 



VE 



^0iir]',x) + 



(57) 

(58) 
(59) 



on the linearized fields. In particular, since under any infinitesimal coordinate 
transformation ()57p the graviton field transforms to: 



for the specific choice 



we obtain: 







1 



dr] a{r]) ^Qo^i[r] ,x) ^— 



^ooiv,^) + H{r]) / dr]' a{r]') ipoo{r]',x) 



VE 



VE i^voii) 



(60) 

(61) 
(62) 

(63) 



Thus, by construction the graviton field (j6HI63p obeys (jl9ll2(J|) . Keep- 
ing in mind the definition the first order temperature fluctuation 
becomes: 



AT, 



■R 



Te 

VR 



e + 



(64) 



K 



dr]' 



Ve 
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Further reduction of uses: 



H{r]) n 



dr]' a{r]') ipoo ,i{r]' , x) 



Ve 



(65) 



together with: 



*o„,o(.,.x-) + !^*„„(,.f) 

a{7]) 

+ H'{r]) r dr]' a{7]') ipM.x) 

■J TIE 



(66) 



A straightforward computation leads to the final form for the first order 
temperature fluctuations: 



AT, 



R , 



T, 



R 



e + 



(67) 



K 

2 



^oo(%,o; 



H 



R 



VR 



dr]' a{r]') ^oo{v\^] 



Ve 



^o'iV- 



■R 



^o'iv 



R 



K 



2 
2 



i'ooiVE , {Vr- VE)e) - K dr]' i^oo^oir]' , {r]R- r]')e) 



Ve 



Vr 



Ve 



dv' -^Ij^oiv' , {VR-V')e) . 



The right hand side of (jU7j) consists of a part associated with the temperature 
fluctuations of the emitting surface plus seven terms. The first two of the 
latter have no angular dependence and belong to the monopole contribution. 
The third term is the dipole contribution while the fourth is the Sachs- Wolfe 
velocity potential term. The spectra of scalar and tensor perturbations that 
are usually reported reside in the fifth (the Sachs- Wolfe potential term) and 
seventh terms respectively. The remaining sixth term is sometimes called the 
integrated Sachs- Wolfe effect. 
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4 The Tensor Power Spectrum 



In 1979, Starobinskii [zo] became the first to calculate the tensor power 
spectrum from what would later be called a model of inflation. Subsequent 
computations were in 1982 made by Rubakov, Sazhim and Veryaskin [21] and 
by Fabbri and Pollock [5^. The definitive result was obtained by Starobin- 
skii in 1985 These calculations all depend upon a normalization for the 
late-time mode functions whose precise determination is our only improve- 
ment. However, we shall also carry out the computation in a slightly different 
fashion. 

The part of (jU7|) relevant to tensor perturbations is: 



AT; 



R , 



R 



Vr 



(68) 



VE 



and can be expressed as a sum over graviton momenta and polarizations: 



AT) 



R , 



R 



J e'eUij{k,s) a{k,s) + (c.c.) } , (69) 



(27r 



where the scalar response function is: 



/i(e, k) 



tE 



UAitj k) j exp 



tR 



ik ■ e dt' a {t') 



(70) 



It is straightforward to compute the expectation value in the presence 
of the state which was empty of gravitons in the distant past: 



a{k,s) \Q) 







(71) 



and obtain: 



n 



ei 



R 



AT, 



R , 



62 



R 



d^k 
(27r)3 



h{ei,k) h*{e2,k) 



X e/ 62 62 



1, 



(72) 



The scalar response function (fTO)) can be explicitly evaluated because the 
physical process occurs entirely during the epoch of matter domination. If 
we assume that the onset of matter domination occurred at a time Im, when 
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the Hubble parameter and scale factor were Hm and qm respectively, then 
at later times: ^ 



Matter 



Hit) 



H, 



M 



1 + i^Hm [t — tM) 



1 + -Hm (t - tM) 



In view of 
Matter 



a{t) = aM 

the dimensionless variable © equals 

3 



x{t,k) = x{tM,k) 



1 + -Hm it - tM] 



(73) 
(74) 

(75) 



In terms of x, the radial component of a lightlike geodesic times k takes the 
form: 

dt' 

Matter =^ k / — — = 2x{tR, k) - 2x{t, k) , (76) 
Jt a[t ) 

and the scalar response function (j7(J|) becomes: ^ 



h{e, k) 



XE 



dx 



d_ 

dx 



UAit,k) 



^2ik-e{xji-x) 



(77) 



Further progress in the evaluation of the scalar response function ()77j] 
requires an explicit form for the mode function. Indeed, the source of our 
improved estimate for the graviton power spectrum is our improved deriva- 
tion of the graviton mode functions ,10j. Since the physical process under 
study involves modes that underwent first horizon crossing at t = ti, the 
relevant form of the mode functions for t > ti is ^T] : 



UA{t,k) 



~zH, r(i-z.) j,,(-|; 



X Ci^(fc) X CiA{k) . (78) 



It consists of three factors, the first of which is the time dependent part: 

/TT J'd'^x) 



r(l - z.) J_.(-|- 



x_\-v 
2q) 



(79) 



2 (2x)i 

2 ' ' 2 ^ ' 

'^During matter domination, the deceleration parameter q{t) is quite well approximated 
by the constant Qm = +\- 

^Henceforth in this section, all quantities refer to the form they take for a matter 
dominated universe. 

''The subscript 1 in a quantity signifies its value at first horizon crossing t — ti. 
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This is a standard result. ^ The remaining two factors in (j78|) represent our 
improvement to the normahzation of the mode functions; Cia depends upon 
the state of the system at first horizon crossing: 



( — -y^ 



(80) 



To get a feehng for Cia, it is important to note that for perfect de Sitter 
inflation it equals one: 



qi = -1 =^ CiAik) = 1 , 
while for a more realistic situation one finds to first order: 
qi{k) = -1 + Aq{k) =^ 



(81) 



CiA{k) 



^il) + ln2 - 1 Aq{k) 



r{z) 



(82) 



The factor depends upon previous evolution. Had there been no evolution 
in q from the initial time ti to ti, its value would be one: 



q{t) = q 



CiA{k) 



(83) 



When - as is the physical case - there is a mild evolution, it results in small 
deviations about (jHBI) whose explicit form is given in the first Appendix. 

In view of ()79j) . we can express ()78j) with its conventional slow-roll nor- 
malization times the two correction factors: 



UA{t.k) 



-iHi 



sin(2a;) cos(2a:; 



Ax' 



With the infrared approximation 
the scalar response function (fTTjl : 



/i(e, k) 



inHi 3 
V2F 72 

sin(2a;) cos(2a;) 

Ax 

Ei(2i(l - 



— iw 



n n 2i w X[ 



X CiAik) X C.Aik) . (84) 
it is possible to exactly evaluate 

sin(2x) cos(2x) 



— w 



4x2 



(85) 



sin(2x) 



— 2i wx 



w )x 



Ax 

- Ei(-2z(l + w) 



X 



xe 



*See, for instance, equation (4.29) of jH]. 
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where, to economize on writing, we have defined w = k ■ e as the cosine of 
the angle between the unit vectors k and e. However, there is no point in 
retaining the full complexity of this result. It is easy to check that the term 
inside the curly brackets falls like for large x. Therefore, potentially 
observable effects must derive from modes which had not yet experienced 
second horizon crossing at the time of emission. This implies -C 1. The 
modes which produce anisotropics within our current horizon volume must 
also have experienced second horizon crossing by the time of reception. Hence 
we can also assume xr':^ 1. It follows that the only significant contribution 
comes from the lower limit, for which we may as well take the limiting form 
relevant to small xe'- ^ 



h{e, k) 



k^ 

1 3 5 
X <! -w 

2 4 



(86) 



-w (1 — w 



, , l+w 
In I I + tiT 

1 — w 



The angular dependence in our expression ()86|) for the scalar response 
function is complicated. However, one can recognize some of the factors as 
spherical harmonics with zenith angle 6 = arccos(e ■ k) and azimuthal angle 
6 = 0: 



1 
2 



4^ 



TT 



Yr 



00 



^1^20 



(87) 



w {1 — w^) 



27r 
105 



32 



It makes sense to decompose the scalar response function into a part depend- 
ing only upon k = \\k\\ and an angular factor B, with the Yqo term in the 
latter bearing unit normalization: 



h{e, k) 



CiA CiA ^ 0(e, k) 



Obviously: 



e(e,fc) 



J2i wxf 



An 



3w' 



- w (I — w'^ 
2 ^ 



In 



l + w 
1 — w 



(89) 



(90) 



^Although our technique has been different, this result seems to agree with Starobin- 
skii's equation (12) |26j . 
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We define the "graviton power spectrum" in terms of tlie radial factor: 



GHl{k) ClAk) \\CUk) 



(91) 
(92) 



Because the hterature abounds with different conventions for this quantity, 
we correspond Vh{k) to the symbol S^ik) used by Mukhanov, Feldman and 
Brandenberger jB], to the variable Vg{k) used by Liddle and Lyth and to 
the quantity A^{k) used by Lidsey et al. [7]: 



Vhik) 



97r 



TT 



V,{k) 



257r 



Al{k) 



(93) 



4 ' 16 ''^ ' 4 

Perhaps the clearest specification of Vh{k) is to state how it enters the tem- 
perature correlation function: 



R , 



AT, 



R , 



R 



62 



R 



°° dk 



Vh{k) 



(fk 
in 



e(ei,fc) Q*{e2,k) 



A C]^ di ^2 ^2 



(94) 



The leading order slow roll result for Vh{k) is typically expressed in terms 
of the value of the scalar potential at horizon crossing. Using (|38|) it can be 
converted to our notation: 

2-gi 



3 



GHf 



(95) 



Our correction factors of 33^) ^lAi^) W^iAik)]]"^ are typically near one 
for slow roll inflation. Note especially the factor ||Cj^(A;)|p, which represents 
the effect of evolution from the beginning of inflation up to horizon crossing, 
as required by the analysis of Wang, Mukhanov and Steinhardt • 

It is elementary to verify that there is no monopole contribution to ()94p 
by fixing one of the two directions, for instance 62, and integrating over the 
other: 



Monopole 



47r ' 



AT. 



AT, 



T 



R . 



R 



T 



62 



R 



(96) 



If we take the 2;-axis to be along the k direction, we can express ei in terms 
of the zenith angle 6 and azimuthal angle 0: 



ei 



sin 6 cos (h , sin 6 sin 6 , cos 6 



(97) 
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The resulting azimuthal integration is simple: 



2-IT 



27r 



(98) 



and the properties of Uij ensure that (j98j) gives a vanishing monopole con- 
tribution (jnni)- 

In a similar fashion, it can be proved that contains no dipole com- 
ponent: 



Dipole 



An 



Seiei in 



AT, 



R , 



AT, 



T 



R , 



R 



T 



62 



R 



, (99) 



where, as for the monopole case, direction 62 has been fixed. 



5 The Scalar Power Spectrum 

The spectrum of scalar perturbations can be computed from the Sachs- Wolfe 
potential term in ^7\i : 



AT 



R , 



sw 



2 ■ipooiVE, {VR-VE)e' 



(100) 



By virtue of (j^^ we have: 



AT 



R , 



T 



R 



SW 



d^k 



UciVE, k) e 



Y{k) + (c.c. 



2 J (27r)3 

In the presence of the state without any scalars in the distant past: 

Y{k)\n) = , 

the temperature correlation function ()30|) becomes: 



(101) 
(102) 



AT. 



AT. 



T 



R 



T 



le2j 



R 



I sw 



i<? f d^k 



4 J (27r)3 



Uc{VE.k)\\' 



X e 



2i {xji-XE)k-{7i -62) 



(103) 



in terms of the dimensionless variable 0- 
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The relevant form of the mode functions is for t > ti jTT] : 

-H, H{t) 



Uc{t, k) 



x< 1 



^2k^ (1 + gi) a{t) 



dt' a{t') [ 1 + q{t') ] 

X CUk) X CUk) , (104) 



where ti, as always, signals first horizon crossing. In analogy with Section 4, 
the normalization factor Cic depends upon the state of the system at ti. It 
is expressed in terms of r{t) - defined in (fH^ - and the parameter qcif)'- 



qc{t) = 



+ 



r{t) 
Hit) 



l + r(t) [l + r{t)y 



(105) 



The expression is: 



V^'^y 2^ qiC 



( 291C ) 



e nc cos(^) ( 1 + ri ) «ic 



lie 



' lie 



qic = qcih) , ri = r{ti) 



(106) 



(107) 



During inflation - in fact, quite generally - the parameter r is typically zero. 
Therefore: 



r{t) = =^ qc{t) = -q{t) =^ Cic{k) = CiA{k) e^i cos(^) 

More generally, if r is small we can write to first order: 

f(t) 



qc{t) = -q{t) + q{t)r{t) + 



Hit) 



(108) 
(109) 



Consequently, as in (jHzj) . we have to first order: 
qi{k) = -1 + Ag(A;) =^ 



Cicik) 



1 + 



Aq{k) + ri 



Hit) 



;iio) 



The other factor, Cic, depends upon evolution from ti to ti. Just like Qa, 
it equals one when q is constant; its general form can be found in the first 
Appendix. 
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Because the physical process takes place entirely during pure matter dom- 



ination: 

Matter 



10 



dt' a{t') [ 1 + q{t') 



a{t)[l + q{t)] 3aE 



Hit) 



2Hf 



■ (111) 



Thus, the mode functions can be expressed as a product of the conventional 
slow roll normalization with the two correction factors: 



Uc{tE,k) 



3Hi 



2 ^2/£3 (1 + q, 

and the temperature correlation function takes the form: 



X Clcik) X Qcik) , (112) 



R , 



R 



le2j 



R 



n 



97rG 



sw 



(27r)3 2A;3 1 + q^ 
x||Cic(A:)f e2^(^«-^^)^-(^i-^2) (113) 



The identity: 



f{k) e 



ik ■ X 



(27r)3 

reduces (jll3p to its final form: 



^ r dkk^ f{k) ^ 



(114) 



ATr ATr 



9G rfA; Hi 



sw 



An Jo k 1 + qi 



C^cik) C,c{k) 



sm 



2(xi? - Xij) II ei - 62 



2(xij - xe) II ei - 62 



■ (115) 



The "scalar power spectrum" is defined by the way it enters the correla- 
tion function between temperature fluctuations observed from directions ei 
and 62: 



ATr,^ ATr 

,ei) ^^(62) 



T, 



R 



n 



/sw -'o /c 



°°dk 



-T-Vswik) —e 



d^k 



{xj{-xe) k-(ei-e2) 



4tt 



(116) 



^°Henceforth in this section, all quantities refer to the form they take for a matter 
dominated universe. 
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Hence, we obtain: 



Vswik) = \\Cwik)r \\C^cik)r . (117) 

We again correspond Vsw{k) to the symbol 6{k) used by Mukhanov, Feldman 
and Brandenberger to the variable Vnik) used by Liddle and Lyth [H], 
and to the quantity A\{k) used by Lidsey et al. [7]: 

^sw(A:) = ^llWf = '-Vn{k) = • (US) 

The leading slow roll result for Vsw{k) is usually expressed in terms of 
the scalar potential and its derivative at the time of horizon crossing. Using 
and (PUI) we can convert this to our notation: 

ii = A ^ ,119) 
4ir 1 + ,1 3(2 - ft + ri)2 ■ ' ' 

Our correction factors of ^^^^z^^^^p^, ||C'ic(^)ll^ HCjclfc)!!^ are typically 
near one for slow roll inflation. Consistent with the analysis of Wang, 
Mukhanov and Steinhardt [TH], there is a factor [[(^^^(fc)!^ which represents 
the effect of evolution from the beginning of inflation up to horizon crossing. 



6 Epilogue 

We have taken advantage of a recent, exact solution for the mode functions of 
scalar-driven cosmology ^1] to re-compute the scalar and tensor power spec- 
tra for anisotropies in the cosmic microwave background. For completeness, 
and to emphasize its inherent gauge invariance, we have also reviewed the 
standard computation of the Sachs- Wolfe effect. The principal new feature 
is our expressions for the normalization factors that were built-up during 
inflation. 

We have not expanded the temperature correlation function in spherical 
harmonics. Nonetheless, since our results take the form of the standard 
normalization times correction factors, it should suffice to simply multiply 
the standard result by these correction factors evaluated at the wavenumber 
appropriate for the l-th multipole moment: 

k = haoHo . (120) 
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The tensor correction factors Cia and are given by (IHU)) and (|132|) respec- 
tively; the analogous scalar factors Cic and Qc by (|1(J6|) and (|133p . 

How observable are the correction factors we have found? Since it is likely 
to require a major effort to detect a non-zero tensor amplitude, the fractional 
improvement we give for this probably does not matter. On the other hand, 
precision measurements of the scalar amplitude might very well be sensitive 
to the structure we provide. The greatest advantage of our formalism is 
not the incremental improvements it offers for the standard, slow roll regime 
but rather its applicability to exotic scenarios that lie beyond the slow roll 
paradigm. We present an example in the second appendix. 

Finally, we disagree slightly with the standard treatment of the tensor 
contribution. The original authors seem to have averaged over graviton po- 
larizations before taking the expectation value. This makes a small but possi- 
bly significant difference in the tensor contribution to the multipole moments 
of the temperature fluctuations correlation function. 
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7 Appendix: The evolution dependent cor- 
rection factors 

A central feature of our exact solutions is the transfer matrix, A4j[t,ti, k). 
There is an / = A transfer matrix for the graviton mode function and an 
I = C one for the scalar mode function. Each of them is the time-ordered 
product of the exponential of a line integral: 



Mi{t,ti,k) = pjexp ^ dt' Ai{t',k) 

= E / / dt2... dt^ Ai{h,k)---Ai{tn,k) . (122) 



The exponent matrix Ai{t, k) vanishes whenever there is no evolution of the 
appropriate qi(t): 

Graviton =^ lAit) = lit) , (123) 

(f) — 

Scalar ^ qMt) = - ! , + ^ — ^ . (124) 

There is a similar dichotomy for the appropriate physical wave number ex- 
pressed in Hubble units, 

k 

Graviton =^ XA{t,k) = , , , , , (125) 
^ ' ^ a{t)H{t) ^ ' 

Scalar xc{t,k) = ^"^^ '^^j . (126) 

With these definitions the exponent matrix takes the form: 



TT 



csc(z/,7r)c.,(-2.) -2zd.,(-^) 



4 \-2tcsc^{iyiTT)K^{-^) -csc(z//7r) c^,(-|^) y 
where the various coefficient functions are: 

^ ^ k rin)Tin-u + l)Tin-2u + l) ' ^'^""^ 



"'^^Recall the definition H^W|I of the parameter r{t). 
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smlz/TT 



1 - In(i^) 



-1)" r(n- i) 



-1 oo 

^ V — 

TT ,^fr^i r(n + z/) r(n + 1) r(n - z/ + 1) ' 



1 



E 



-1 " r(n + u- i 2 



1\ ^2n+2;/ 



n + u) 



(129) 
(130) 

(131) 



2v^ ^0 r(n + 2z/) r(n + z/ + l) r(n + l) 
and we have defined: 

We can now give precise definitions for the evolution dependent normal- 
ization factors: 

CiA{k) = M"^ {h, ti, k) + M'^{ti,ti,k)e'fi sec{f,) , (132) 
Cic(k) = M^(^ {h,ti, k) + M^J{h,U,k) e-'^ sec{^) . (133) 

The subscript i denotes the initial value of the respective parameter. Since 
during inflation one typically has: 



(134) 



it ought to be a very good approximation to simply take the first several terms 
of the series expansion of the transfer matrix in estimating these corrections: 



MY 



1 + 



dt-fi{t) + / dt dt'\-fi{t)-fi{t')-Si{t)f3i{t') 



Mf ~ -i r dt 5i{t) - i r dt f dt'[^i{t) 5i{t') - 5i{t) ^i{t') 
Mf - -i r dtf3i{t) - i r dt f dt'[f3i{t)^i{t')-^i{t)f3i{t') 

J ti J ti --fti 

dt ji{t) + r dt t dt' [ 7^(0 -fi{t') - Pi{t) 6i{t') 



Mf ~ 1 - 

where the coefficient functions are: 

6i{t,k) = 



2 sin (z^/vr) 
4sin(z/77r) 



duA-f,) ■ 



2 -^''I V qi > 



(135) 
(136) 
(137) 
(138) 

(139) 

(140) 
(141) 
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8 Appendix: Ultra Slow Roll Inflation 



Consider an inflaton potential like that depicted in Fig. 4 and suppose infla- 
tion begins with the scalar to the right of the flat portion. Once the scalar 
rolls into the flat region its background equation of motion becomes: 

^o + 3Hipo = . (142) 

This can be integrated to give an exact expression for the scalar's time deriva- 
tive in terms of its value at the beginning of the flat region: 

If the scalar has enough kinetic energy it can roll through the flat region, 
and then on down its potential. The condition for this to happen is: 

We shall assume this and study the scalar power spectrum for modes which 
experience first horizon crossing while the scalar is on the flat section. 

In the flat region all derivatives of the potential vanish, so all of the slow 
roll parameters are zero. Although the scalar is rolling ever more slowly 
- hence the name - this is a situation in which the conventional slow roll 
approximation completely breaks down. In fact the slow roll prediction (jll9|) 
for the scalar power spectrum actually diverges! The difficulty of reconciling 
this with a system which is approaching a pure de Sitter phase was the 
occasion of much reflection by Grishchuk |27j. We shall see that Vsw{k) is 
finite, but that it can become quite large. 

By adding the background Einstein equations ()35II36|) and then substi- 
tuting (|143p one finds: 

l+,W.4.G||^4.G(|Ly(^)\ ,145) 

During inflation the deceleration parameter is typically near —1, but the fact 
that it approaches this value exponentially fast during the ultra slow roll 
phase makes a crucial change in the parameter r{t) defined in (jH^ : 

r{t) ^ J^JtH^/^^) = -^-§-2 = + ■ (146) 
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Figure 4: The scalar potential associated with a phase of ultra slow roll inflation. 

In the region y^e < < V?/ the potential is exactly flat with value Vf. 

Although r(t) is near zero for typical models of inflation, we see that it is 
nearly —3 during the ultra slow roll phase. It is simple enough to obtain an 
exact expression as well for its derivative during this phase: 

Note that this quantity is nearly zero, both for typical inflation and during 
ultra slow roll phase. 

It is now straightforward to evaluate our factor Cic{k) that depends upon 
the system's state at horizon crossing. Substituting in ()105|) gives the follow- 
ing result during the ultra slow roll phase: 

n (f) - ^ 2(2-g)(l + g) _ 1 A 

Although qc{t) ~ +1 in typical models of inflation, we see that it rapidly 
approaches — | during the ultra slow roll phase. Evaluating (jl(J6|) for qic = 
— I and ri = —3 gives: 



qic -- 



1 „/ 3x 1 



40? V 2^ 3 

ri = — 3 V 



r -o = o • (149) 
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To estimate the evolution-dependent factor Cic{k) we make the reasonable 
assumption that the system goes suddenly from qc '^^ +1 to qc — —\- In 
this case the transfer matrix is determined by matching the mode functions 
and their first time derivatives at the onset of the flat region: 



(-iJ^^_{-xj), J,_{-X}))(-^g -2^22) = (iJ^i-Xf), J_^_{-Xf)), (150) 
(-./,(-x^),4(-x^))(t;ji = {^A{-XJ),J\{-Xf)). (151) 

The matrix elements needed for the scalar power spectrum are: 



c 



Ml' 



m Xf 

2 

nxf 



-J_d-Xf) JU-Xf) + J'_^{-Xf) Ji{-Xf) 
-J_5(—Xf) J'_i{—Xf) + J'_5{—Xf) J_i{—Xf) 



Substituting 


in ()133|) with qic 


= +1 we 


Qcik) = 


Ml' + M'J , 











2 



-Xf) 



lJ's{-Xf) 



3 3i 
■1 + — + ^ 



X, 



— sm(x/j H — 3 cos(x/J 



X 



f 



X 



(152) 
(153) 

(154) 
(155) 

(156) 
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Because first horizon crossing occurs after the scalar has rolled onto the 
fiat region we can assume Xf > 1. It is not safe to assume x 1 because 
some modes will experience horizon crossing soon after the ultra slow roll 
phase begins. The power spectrum of these modes will deviate much more 
from scale invariance than is typically the case. Although the fiat region must 
be narrow enough that the scalar can roll across, this process can be tuned to 
require an arbitrarily long time. For modes which experience horizon crossing 
long after the onset of the ultra slow roll phase, one can assume x/ ^ 1, in 
which case: 



C^c{k) 



(157) 



X;>1 =^ 

We constructed this model as an exotic system in which the slow roll 
paradigm completely breaks down. However, it has two other properties 
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In accordance with the definition © , a;/ 



Hf 
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worthy of note. The first is that, although our prediction (jll7|) for the scalar 
power spectrum remains finite, it can become quite large owing to the inverse 
factor of (1 + qi). We have seen from ()145p that (1 + q{t)) approaches zero 
exponentially fast. It seems inevitable that back-reaction must eventually 
become significant if the ultra slow roll phase is protracted. 

The second interesting property of this model is that the anisotropics 
generated during the ultra slow roll phase are entirely due to scalar kinetic 
energy. The potential is completely fiat so the only possible fiuctuations 
derive from the gravitational response to kinetic energy. This is usually 
dismissed as negligible but we have just seen that it can drive an enormously 
strong effect as the system approaches de Sitter infiation. This suggests that 
one might expect a similarly strong effect from gravitons - the combination 
of two of which can produce a scalar - if the computation were carried to 
next order in the weak field expansion. 



30 



